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- Finally, similar curves were drawn for System T—mean 
photographic R. A., and the results were decidedly more consistent 
than for C. 

23. From these investigations, which have been tested and 
varied in every way that can be thought of, the following results 
■emerge:— 

System T has a fairly uniform magnitude-equation of about 
+ o s ‘oi5 per magnitude throughout list I. and the first half of 
list II. ; for the second half there is little material. 

System L and the independent adopted systems have a 
magnitude-equation of about +o s- 02o for list I. and +o 8, oi5 
for list II. 

System C has a magnitude-equation of about —o s- oi7 for 
about half list I. and of — o s *oo2 for the remainder and for 
list II. 

It seems to me, therefore, that so far as magnitude-equation 
is concerned System T is little better than System L, while 
■System C is worse, because it changes so suddenly. From this 
point of view there would be little advantage in adopting either 
T or C as the standard system for the reduction of all the 
photographs. 

My sincere acknowledgments are due to Miss Julia Bell, who 
has carried out very skilfully the numerical work summarised in 
this paper. A large portion of the expense has been borne by a 
grant from the Royal Society Government Grant Fund. 

Cambridge Observatory: 

1906 June 7. 


On some Points connected with the Determination of Orbits. 

By H. C. Plummer, M.A. 

1. The difference between the circular measure and the sine 
of an angle is an expression which occurs in several of the most 
important formulae relating to elliptic motion. When the angle 
is large there is no difficulty in calculating this difference accu¬ 
rately, but when it is small ordinary logarithmic tables will not 
give the required accuracy without the help of some special 
device. Hence auxiliary tables have been published in a variety 
of forms, according to the purpose for which they are designed. 
It is possible to restrict the compass of such tables by making 
them merely supplementary to the ordinary logarithms by the 
use of some artifice. Of this nature is Tietjen’s * formula, 
which may be expressed in the form 

£—sin £=-§- B sin 3 i £ (sec -Je) 2 ' 4 ... ... (1) 

* A.N. 1463 ; also Watson’s Theoretical Astronomy, p. 343. 
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for which a very small table, giving log B as far as e=6o°, is 
necessary. But the use of any special table can be avoided. 
Thus I have pointed out elsewhere * that the function 

i e 3 / (sec ^ £ )*4-4 .(2) 

is a remarkably close approximation to the value of t—sine; 
indeed, over the range of e from o° to 70° the difference between 
the logarithms of the two expressions does not exceed 2 in the 
seventh plaoe. A greater error may be caused by the occurrence 
of the power 14*4. It may be noted that if e is less than 36° 
the logarithm of the required secant can be obtained from 
Schron’s tables (as the difference between S and T) to eight 
places. Within this range at least the approximation may be 
considered practically perfect. 

2. A function of such frequent occurrence that it has been 
made the subject of special tables is 

Q( £ ) = l( f — sin £ )/sin3£s. 

This can be expressed by Tietjen’s formula thus : 

Q(«) = B(secif-) 2 - 4 . (3) 

or by using the approximate expression (2) 

Q( £ ) = (4 £ / sin § £ ) 3 ( sec ts £ ) -14 * 4 . ( 4 ) 

provided £<70°. How the nature of the function Q(e), which 
is unity when £=o°, suggests that it may be represented very 
approximately by the form (sec be) a , a and b being two constants 
so chosen as to give the best possible result. For practical 
purposes, however, b must be a very simple factor. This limita¬ 
tion precludes a quite satisfactory approximation: no great 
improvement can be made on the trigonometrical factor in (3). 
But the case is different if we seek to represent Q not by a 
single term (sec be) a , but by a product of two such terms. This 
is shown by considering the logarithmic expansions of the 
functions involved. 

3. In the first place, if p is the modulus of common loga¬ 
rithms, 

log (e —sin e)—3 log £ 4- log 6 

I68OO £4 4 " >i6 I 000 £ ^d" 3 I 04640000 £ ^ • ‘ •] (5} 

We have also 

2 2 * 

log sin X — log x— uS —. x 2i . (6) 

2 % | 2 % 

log sec x — /iS . (7) 

21 2% 

** ' 


* A.N 3883 
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where B„ B 3 , . . . are Bernoulli’s numbers, of which the first 
four are £, and By means of (5) and (6) we find 

log Q(£)=p [^£ 2 + j i2o~o f ^ “ 1 " i 60 00 t6 4 " 206 of60 ooo £8 + •* •] ( 8 ) 

The object is now to combine two series of the form (7) in such 
a way as to reproduce the first terms of (8). Theoretically an 
attempt might be made to obtain the first four terms; but if 
simple multiples of e only are admitted, a less perfect agreement 
must suffice. It is indeed remarkable enough that with the 
multiples £ and | of e three terms can be made to coincide. 
Thus the application of (7) gives 

log sec ie— log sec 

” "h I I^OO^d - 5 6000 £ ^”l"4fr3 84000^8 "!"••• 

= log Q ( £ ) 455 dgCOOOcA 6 “ • * * 

The outstanding part of the fourth term amounts to 1 in the 
seventh place when £=82° and the succeeding terms are still 
smaller. Practically throughout the first quadrant we may 
write 

log Q (e) = log sec ^ ^ log sec £e 

= [0*5454132] log sec i£—[0-3978408] log sec |-£ ... (9) 

The range is greater than is actually necessary, since Q can 
be calculated directly when £ has a fairly large value. 

4. As an example of the use of such formulae let us consider 
the method of Gauss for determining an orbit when two helio¬ 
centric positions are known. This involves the solution of the 
equations 

y 2 = m/ (l+ sin 2 £ g) .(10) 

y 3 — y 2 — m( 2 g —sin 2g) j sin3<7 = £mQ,(2g) ... ... (11) 

where y is the ratio of the sector to the triangle, 2 g is the 
difference of the eccentric anomalies, and l and m are given 
quantities. These equations can be solved by trial without 
using special tables such as have been given by Gauss. The 
natural course is to choose an approximate value of 2 g (in the 
absence of more precise knowledge, 2 f, the difference of the true 
anomalies, may be taken) and to deduce y by means of (n). 
Then (10) will give an improved value of g, with which the pro¬ 
cess can be repeated. 

Let y — 1 /c sinha. Then, by (11), 

4mQc3 sinh 3 a + 30 sinha—3 = o, 

which can be compared with 

4 sinh 3 a + 3 sinha—sinh3a = o, 


] 
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giving mQc 2 = i, c sinl^a = 3. 

Now if e 3 “ = cot ±ft, e d =cot -|y 

sinh3a = cot / 3 , sinha = cot y. 

For (10) we may write 

sin 2 \g =■ l tan 2 0, y 2 = m cos 2 c> / 7 . 


Then the calculation from an assumed value of g to a better 
approximation is reduced to the following system of equations :— 


with 


cot /3 = 3m*Q* > 

tan 3 iy = tan -|/3 
cos h = tan y - 
sin -|gr = tan 3 
y = cos <)> 


... (12) 


the last step being made only after satisfactory values of g and d 
are known. 

5. The application of the method will be made clearer by a 
numerical example. The following * is chosen : t'—t, the interval 
of time, = 100 days, logr = 0*221 6050, logr'= 0*209 9050, 2/, 
the angle between the two radii, = 44°25 , 48 // *oo ; whence 

m £) 2 /8 cosy(rr')= = [9*021 2961] 

l =(r+r')/4 cos/(rr') } — \ = [8*603 5 66 3] 


The complete calculation by means of Bremiker’s 6-figure 
logarithms is given below (see p. 495). 

Now if, as in the given calculation, one assumed value x T 
leads to a closer approximation x 2 , and x 2 similarly leads to x 3 , 
the approximate correction to x 3 is 

X~~~X 3 — (^s* - ’^j 2 / {(^2*”*^ 2 )} ••• (^3) 

w hich, applied to the successive values of •§<?, gives 
ig—iffz = (63" , 8) 2 ^-3658' / = + i"-n. 

The calculation of y can now be completed thus :— 


hr 

0 in 

12 5 181 

8 

O / // 

46 x6 6*2 

L. sin \g 

9*321018 

L. cos 5 

9^39654 

£L .1 

9-301783 

i L .m 

9-510648 



-*L .1 

0-698217 

L.tan 5 

0*019235 

L .y 

0-048519 


* For the same example worked by Gauss’ special tables see Bauschinger’s 
Tafeln tur Theoretischen Astronomic , p. 26. 
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Calculation referred to on p. 494. 


2 9 

O / // 

44 25 48 0 

0 / it 

48 25 23-2 

0 

O / U 

45 9 20*0 

0 / // 

45 156-2 

i .9 

II 6 27*) 

12 6 20-8 

¥ 

22 34 40*0 

22 30 58-1 


7 24 18 0 

8 4 139 

L. tan |(8 

9*618890 

9-617571 

L. sec \g 

0-008213 

0*009766 

L. tan 

9*872963 

9-872524 

L. sec \g 

0-003637 

O*004323 

47 

0 / // 

36 44136 

0 / a 

36 42 33 6 

LL. sec \ 

7-560743 

7-635785 

y 

73 28 27-2 

73 25 7-2 

const. 

0-397841 

0-397841 




L.(i) 

7 958584 

8-033626 

L. tan 7 

0-527678 

0*526136 




4L .1 

9-301783 

9-301783 

LL. sec \g 

7-914502 

7-989717 

5 L. Q 

0-009872 

0-011741 

const. 

o-5454»3 

0-5454I3 




L. (2) 

8‘4599I5 

8-535I30 

L.cos 8 

9-839333 

8839660 

( 2 ) 

0-028835 

0-034287 

8 

0 / n 

4618 32-3 

0 It 

46 16 3-6 

(0 

0-009090 

0-010805 

L.tan 8 

0-019851 

0-019224 

L.Q 

0-019745 

0-023482 

4L.J 

9-301783 

9-301783 

L- 9 

0-954243 

0-954243 




L. m 

9*02I296 

9-021296 

L. sin %g 

9-321634 

9-321007 

sum 

9-995284 

9-999021 

4 9 

0 i n 

12 6 20*8 

Oil/ 

12 5 170 

L. cot J8 

9-997642 

9-999511 





The value found by 7-figure logarithms and Gauss’ tables is 
log y = 0 0485191. The accuracy of the result is naturally 
greater than can be expected in general from a 6-figure com¬ 
putation* but the above example is merely illustrative of the 
method, and a more accurate calculation is therefore unnecessary. 

6. It may be remarked that a still rougher computation may 
be of use as a preliminary step. Thus 4-figure logarithms are 
enough to give \g in this case with an error not greater than 1'. 
But for the effective application of (13) it is necessary to make 
the two successive calculations with nearly the same order of 
accuracy. This formula, which is quite general, can easily be 
proved on the assumption that only the first powers of the errors 
committed are of importance. The general problem is to solve 
the equations 

y = p(x), y = q(x) s 

where p and q are functions which involve algebraic or transcen¬ 
dental forms, by assuming a value x z and deducing x 2 , which in 
turn is made a starting-point to obtain x 3 . The two equations 
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may be represented graphically by the curves P,;CP 2 and 
Q 2 CQ 3 , the coordinates of whose intersection C are to be found. 
The method may then be translated thus : (1) with the assumed 
abscissa x x find the ordinate y x of P„ (2) with this ordinate find 
the abscissa x 2 of Q 2 on the second curve, (3) with the abscissa x 2 
find the ordinate y 2 of P 2 on the first curve, (4) with the ordi¬ 
nate y 2 find the abscissa x 3 of Q 3 on the second curve. It is 
then clear that the meaning of the correction given by (13) is to 
give the point R in which the straight lines P Z P 3 and Q 2 Q 3 
intersect. This graphical representation shows clearly the 
general nature of the process, and suggests in an interesting 
way the conditions on which its effectiveness depends. These 
conditions involve the slope and curvature of the lines, and are 
well illustrated by the special problem discussed above. 



where 


sin 1 


k(t' — t) = («—sin e) — (S — sin S ) 

r +/+ c , rin ^r +r'- o 
4 a 4 a 


... (14) 

••• (15) 


I have shown * that when the eccentricity is moderate and 
the interval between the observations is not unreasonably large, 
this system of equations is capable of a convenient solution and 
gives a practical method of determining the elements of an orbit. 
The usefulness of the theorem has been more generally recog¬ 
nised in cases of large (nearly parabolic) eccentricity. Now it 
is obvious that a direct calculation of t'—t will not be accurate 
when the mean distance a is large and consequently e and 2 are 
small, especially if the chord c is also small in comparison with 
the radii. Hence Marth f has derived elaborate expansions in 
series and calculated extended tables from the results. The 
value of these tables is well known. Yet, though they are 

* M . N . vol. lxiii. p. 147. f A . N . 1557-1560. 
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doubtless convenient, they are not indispensable. It seems 
•worthy of notice that complete accuracy can be attained by the 
ordinary forms of calculation. 

8. With any assumed value of a, e and 8 can be calculated 
by (15); moreover, the difference of these equations gives 

sin|-(6—8) sin-§■ (e + 8) = 0/2®. ... (16) 

whence the difference e —8 can also be calculated with all neces¬ 
sary accuracy. Now since 

2sini(£ — 8)—(sinf—sin 3 ) = 2sini(e—8) {1 —cosf (e + 8) } 

= c/a. tani(e + 5 ) 

the equation (14) can now be written 

k(t'—t) — a*c tan ^(e + 8) 

+ 2 a t (A(£- 5 )-sin£ (e-8)} ... (17) 

We now have a sum of two positive terms, each of which can be 
calculated without difficulty. By (16) this becomes 

k{t' — t) = a^c tan -J(e + 3 ) + sin ~ 3 ^(c + 8). A ... (18) 

where 

A = 6 .{:§■(£—8)—sin-§..(e-^ 8 )}/.sin 3 |^f—8) 

and tends to the value 1 when £—8 is small. It is always 
possible to calculate log A with the help of (2) or (9). Frequently 
the second term in (18) is small compared with the first term, 
and A is not required with great accuracy,. An approximation 
which may be useful can be obtained by the method of § 3. 
This is 

log A = 3'6 log sec ^ («—S) . (19) 

which exceeds the true value by 1 in the seventh place when 
e—8=12°, and will therefore in general be good enough in the 
case of an ellipse which is nearly parabolic. 

9. It is of interest to notice the corresponding form of Euler’s 
equation for parabolic motion. We have a infinite and 8 =£ = o, 
but at the same time the finite limits 

ae 2 = r+r' + c, aci z = r+r'—c. 

Since now A = 1, (17) becomes 

k(t '—£) = ■JcR+-|c 3 /R 3 ... ... ... (20) 

where 

R = (r + r 1 -j-c) i + (r+r / — c)K 

This, unlike the ordinary form, is suitable for the calculation of 
t’—t however small c may be, although not so convenient as the 
trigonometrical transformation employed by Encke.* 

* Bcrl . Jahrbuch , 1833, p. 268 
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10. The numerical application of (18) is in practice very 
simple and convenient. For example, let r = 1*5, r' = 1*51, and 
c = 0*15, and let us calculate the time-interval corresponding to 
these values * of a : 1-55, 10, 400 and infinity. We obtain the 
following results :— 


a. 

Period. 

Years. 

h • 

0 / // 

15 . 

0 /n 

42 46 47-55 

v-u 

Days. 

i *55 

i -93 

45 3316-26 

10-549300 

10 

316 

1619 26-55 

IS 30 33 49 

'7-865279 

400 

8000 

2 32 49-630 

2 25 23-236 

7-570711 

00 

CP 

... 

... 

7-563420 


In all these cases the second term on the right of (18), or (20) 
in the case of the parabola, could be calculated by 4-figure 
logarithms. Further, the value of log A was 0-0005 in the first 
case, and negligible for the other two elliptic orbits. The 
form (18) will always give accurate results, and the advantage 
which it offers becomes most conspicuous in those cases which 
have been considered to present the greatest difficulty. 

University Observatory, Oxford : 

1906 June 4. 


Errors of Tabular Place of Jupiter , from, Photographs taken with 
the Astrographic 13-inch Refractor of the Royal Observatory, 
Greenwich. 

(Communicated by the Astronomer-Royal.) 

In communicating the results of measures of photographs of 
the sixth and seventh satellites of Jupiter it was pointed out 
that in deducing the position angle and distance the error of the 
tabular place of Jupiter had been neglected. 

To eliminate this error, and also that arising from any 
systematic error of the catalogues employed, a series of photo¬ 
graphs of Jupiter was taken with the astrographic 13-inch 
refractor, with exposures only just long enough to give good 
measurable images of the known stars, i.e. 30 s . 

In all ten photographs, taken between 1905 November 3 
and 1906 February 15, have been selected and measured. Four 
images each of Jupiter and of about twelve stars were measured 
on each plate, the positions of the stars being derived from the 
Astronomische Gesellschaft Catalogue (Berlin zones). The deduced 
positions of Jupiter are thus affected by any systematic error of 
the catalogue and the error due to twenty-five years’ unknown 
proper motion of the stars ; but, as the positions of the satellites, 
deduced in the same manner, are affected by the same errors, it 

* For two methods of calculating t' — t for the case a —10 see Bauschinger’s 
Tafeln, pp. 32, 33. 
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